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Abstract 

In this paper, collaborative use of relays to form a beamforming system and provide physical-layer security is 

<N 

investigated. In particular, decode-and-forward (DF) and amplify-and-forward (AF) relay beamforming designs under 

M 

total and individual relay power constraints are studied with the goal of maximizing the secrecy rates when perfect 
channel state information (CSI) is available. In the DF scheme, the total power constraint leads to a closed-form 
solution, and in this case, the optimal beamforming structure is identified in the low and high signal-to-noise ratio 
(SNR) regimes. The beamforming design under individual relay power constraints is formulated as an optimization 

C/3 ' 

O . problem which is shown to be easily solved using two different approaches, namely semidefinite programming and 

second-order cone programming. A simplified and suboptimal technique which reduces the computation complexity 

> 

VjQ ' under individual power constraints is also presented. In the AF scheme, not having analytical solutions for the 

00 ' 

■ optimal beamforming design under both total and individual power constraints, an iterative algorithm is proposed to 

, numerically obtain the optimal beamforming structure and maximize the secrecy rates. Finally, robust beamforming 

o 

designs in the presence of imperfect CSI are investigated for DF-based relay beamforming, and optimization 
frameworks are provided. 

Index Terms: amplify-and-forward relaying, decode-and-forward relaying, physical-layer security, relay beam- 
forming, robust beamforming, second-order cone programming, secrecy rates, semidefinite programming. 



I. INTRODUCTION 

The broadcast nature of wireless transmissions allows for the signals to be received by all users within 
the communication range, making wireless communications vulnerable to eavesdropping. The problem of 
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secure transmission in the presence of an eavesdropper was first studied from an information-theoretic per- 
spective in HI where Wyner considered a wiretap channel model. Wyner showed that secure communication 
is possible without sharing a secret key if the eavesdropper's channel is a degraded version of the main 
channel, and identified the rate-equivocation region and established the secrecy capacity of the degraded 
discrete memoryless wiretap channel. The secrecy capacity is defined as the maximum achievable rate from 
the transmitter to the legitimate receiver, which can be attained while keeping the eavesdropper completely 
ignorant of the transmitted messages. Later, Wyner's result was extended to the Gaussian channel in [3] and 
recently to fading channels in H| and [[5]|. In addition to the single antenna case, secrecy in multi-antenna 
models is addressed in |[6l and 0. One particular result in [6] and [7] that is related to our study is that 
for the MISO secrecy channel, the optimal transmitting strategy is beamforming based on the generalized 
eigenvector of two matrices that depend on the channel coefficients. Regarding multiuser models, Liu et al. 
[8] presented inner and outer bounds on secrecy capacity regions for broadcast and interference channels. 
The secrecy capacity of the multi-antenna broadcast channel is obtained in 0. 

Having multiple antennas at the transmitter and receiver has multitude of benefits in terms of increasing 
the performance, and provides the potential to improve the physical-layer security as well. Additionally, 
it is well known that even if they are equipped with single- antennas individually, users can cooperate to 
form a distributed multi-antenna system by performing relaying lfT3~l - lfT5l . When channel side information 
(CSI) is exploited, relay nodes can collaboratively work similarly as in a MIMO system to build a virtual 
beam towards the receiver. Relay beamforming research has attracted much interest recently (see e.g., 
lfT6l - ll2"0l and references therein). The optimal power allocation at the relays has been addressed in ifTTl 
and lfT8l when instantaneous CSI is known. In lfT9l . the problem of distributed beamforming in a relay 
network is considered with the availability of second-order statistics of CSI. Most recently, Zheng et al. Il20l 
have addressed the robust collaborative relay beamforming design by optimizing the weights of amplify- 
and-forward (AF) relays. They maximize the worst-case signal-to-noise ratio (SNR) assuming that CSI 
is imperfect but bounded. Transmit beamforming and receive beamforming strategies have been studied 
extensively for over a decade. A recent tutorial paper 021 provides an overview of advanced convex 
optimization approaches to both transmit, receive and network beamforming problems, and includes a 
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comprehensive list of references in this area. 

Cooperative relaying under secrecy constraints was also recently studied in |[2T | - ll23l . In EH, a decode- 
and-forward (DF) based cooperative protocol is considered, and a beamforming system is designed for 
secrecy capacity maximization or transmit power minimization. For amplify-and-forward (AF), suboptimal 
closed-form solutions that optimize bounds on secrecy capacity are proposed in [22]. However, in those 
studies, the analysis is conducted only under total relay power constraints and perfect CSI assumption. In 
this paper, we investigate the collaborative relay beamforming under secrecy constraints in the presence 
of both total and individual power constraints with the assumptions of perfect and imperfect channel 
knowledge. 

More specifically, our contributions in this paper are as follows: 

1) In DF, under total power constraints, we analytically determine the beamforming structure in the 
high- and low-SNR regimes. 

2) In DF, under individual power constraints, not having analytical solutions available, we provide an 
optimization framework to obtain the optimal beamforming that maximizes the secrecy rate. We use 
the semidefinite relaxation (SDR) approach to approximate the problem as a convex semidefinite 
programming (SDP) problem which can be solved efficiently. We also provide an alternative method 
by formatting the original optimization problem as a convex second-order cone programming (SOCP) 
problem that can be efficiently solved by interior point methods. Also, we describe a simplified 
suboptimal beamformer design under individual power constraints. 

3) In AF, we first obtain an expression for the achievable secrecy rate, and then we show that the optimal 
beamforming solution that maximizes the secrecy rate can be obtained by semidefinite programming 
with a two dimensional search for both total and individual power constraints. 

4) Two robust beamforming design methods for DF relaying are described in the case of imperfect CSI. 
The organization of the rest of the paper is as follows. In Section HH we describe the channel model and 

study the beamforming design for DF relaying under secrecy constraints. Beamforming for AF relaying 
is investigated in Section |inj In Section [TV] robust beamforming design in the case of imperfect CSI is 
studied. Numerical results for the performance of different beamforming schemes are provided in Section 
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Fig. 1. Channel Model 



|V] Finally, we conclude in Section |VIJ 

II. Decode-and-Forward Relaying 

We consider a communication channel with a source S, a destination D, an eavesdropper E, and M 
relays {R m }^ =1 as depicted in Figured] In this model, the source S tries to transmit confidential messages 
to destination D with the help of the relays while keeping the eavesdropper E ignorant of the information. 
We assume that there is no direct link between S and D, and S and E, Hence, initially messages transmitted 
by the source are received only by the relays. Subsequently, relays work synchronously and multiply the 
signals with complex weights {w m } and produce a virtual beam point to the destination. We denote the 
channel coefficient between the source S and the m th relay R m as g m 6 C, the channel coefficient between 
R m and the destination D as h m G C, and the channel coefficient between R m and eavesdropper E as 

It is obvious that our channel is a two-hop relay network. In the first hop, the source S transmits x s to 
the relays with power i?[|x s | 2 ] = P s . The received signal at R m is given by 

where rj m is the background noise that has a complex, circularly symmetric Gaussian distribution with 
zero mean and variance of N m . 
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In the second hop, we employ decode-and-forward transmission scheme. In this scheme, each relay 
first decodes the message x s and normalizes it as x' s = x s /\/~P~ s . Subsequently, the normalized message is 
multiplied by the weight factor w m by the m th relay to generate the transmitted signal The 
output power of the m th relay R m is given by 



E[\x r \ 2 } = E[\w m x'f} = \w m \ 2 . (2) 

The received signals at the destination D and eavesdropper E are the superpositions of the signals 
transmitted from the relays. These signals can be expressed, respectively, as 

M 

Vd = ^2 h rnW m x' s + n = h f wx' s + n , and (3) 

m=l 
M 

y e = ^2 z mW m x' s + m = z W s + n x (4) 

m=l 

where no and n\ are the Gaussian background noise components at D and E, respectively, with zero mean 
and variance N . Additionally, we have defined h = [h*, ....h* M ] T ,z = [z\, ....z* M ] T , and w = [wi, ...wm] t 
where superscript * denotes conjugate operation, and (-) T and (-)t denote the transpose and conjugate 
transpose, respectively, of a matrix or vector. The metrics of interest are the received SNR levels at D and 
E, which are given, respectively, by 

r = 12^=i n m w m \ ^ = \^ m=1 z m w m \ 
N N Q 

It is well-known that given the channel coefficients, the secrecy rate R s over the channel between the 
relays and destination is (see e.g., 0) 

R s = I{x' s ] y d ) - I{x' S i Ve) (6) 

= log(l + r rf )-log(l + r e ) (7) 



log - """" (8) 

=1 z m w ri 



where /(•; ■) denotes the mutual information, and x' s is Gaussian distributed with zero-mean and E^a/J 2 ] = 
1. Coding strategies that achieve the secrecy rates involve randomization at the encoder to introduce 



uncertainty to the eavesdropper. Secrecy coding techniques are discussed in detail in Q]| - Q. Practical 
coding schemes for secure communications have been studied in ifTOll and [fTTT] for certain special cases of 
the wiretap channel. It is important to note that we assume in the decode-and-forward scenario that the 
relays use the same secrecy codebook and transmit the same signal x' s simultaneously. We further note that 
we throughout the text are interested in beamforming vectors that satisfy for given channel coefficients the 
inequality, N + \ J2 m =i ^rnW m \ 2 > N + \ J2m=i z mW m \ 2 - If there are no such beamforming vectors and 
the ratio inside the logarithm in ([8]) is less than 1 , then the secrecy rate, by definition, is zero meaning that 
secure transmission cannot be established. The beamforming vectors which lead to zero secrecy capacity 
are not of interest. 

In this section, we address the joint optimization of {w m } and hence identify the optimum collaborative 
relay beamforming (CRB) direction that maximizes the secrecy rate given in d8]). Initially, we assume 
that the perfect knowledge of the channel coefficients is available. Later, in Section [TV] we address the 
case in which the channel coefficients are only imperfectly known. We would like to also remark that the 
secrecy rate expression in © in a fading environment represents the instantaneous secrecy rate for given 
instantaneous values of the channel fading coefficients. Hence, in such a case, our formulation considers 
the optimization of {w m } in order to maximize the instantaneous secrecy rates. 

A. Optimal Beamforming under Total Power Constraints 

In this section, we consider a total relay power constraint in the following form: ||w|| 2 = < P T . 
The optimization problem can now be formulated as follows: 




(9) 



wt(^I + hht)w 



(10) 



wtw<P T wt(^I + ZZT)w 

wt(iV I + P T hht)w 

log max — — 77 — 

w t w <p T wT(jVoI + PtZZ'JW 

log A max (iVoI + P T hh\ N Q I + P T zJ) 



(11) 



(12) 
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where A max (A, B) is the largest generalized eigenvalue of the matrix pair (A, B) Lj. Hence, the maximum 
secrecy rate in (fT2|) is achieved by the optimal beamforming vector 

w opt = qu (13) 

where u is the eigenvector that corresponds to \ maiX (N I + P T hh^ , N I + P T zz^) and <? is chosen to ensure 
w lpt w o P t = Pt- Note that in the first-hop of the channel model, the maximum rate we can achieve is 

R x = min log (l + \EAf±) . (14) 



m=l,...,M y jV,.. 

Since we want all relays to successfully decode the signal transmitted from the source in the DF scenario, 
the rate expression in (fl?)) is equal to the minimum of the rates required for reliable decoding at the relays. 
Hence, the first-hop rate is dictated by the worst channel among the channels between the source and the 
relays. 

The overall secrecy rate is 

Rdof, a = mia(R 1 ,R s ). (15) 

Above, we observe that having a severely weak source-relay channel can significantly degrade the perfor- 
mance. In these cases, other forwarding techniques (e.g., amplify-and-forward) can be preferred. Throughout 
the analysis of the DF scenario, we will not explicitly address these considerations and we will concentrate 
on the secure communication between the relays and the destination. Hence, we will have the implicit 
assumption that the the source-relay links do not constitute a bottleneck for communication. 

Next, we provide some remarks on the performance of collaborative relay beamforming in the high- 
and low-SNR regimes. Optimal beamforming under total power constraints is studied in detail in |f2Tfl and 
|[22|. However, these studies have not identified the beamforming structure at low and high SNR levels. 
For simplicity, we assume in the following that the noise variances at the destination and eavesdropper are 
N = l. 

'For a Hermitian matrix A € c ix " an( j positive definite matrix B G C nxi \ (A, iji) is referred to as a generalized eigenvalue - eigenvector 
pair of (A, B) if (A, ip) satisfy Axjj = XBip (28]. 



1) High-SNR Regime: In the high SNR scenario, where both P s , P T — >■ oo , we can easily see that 



Km - logP s ) = min log(|5f m | 2 /iV, 



Ps— >00 



m=l,...,M 



(16) 



From the Corollary 4 in Chapter 4 of 0, we can see that 



lim (R s — \og(P T )) = log(max Ih^j 2 ) 



(17) 



where xp is a unit vector on the null space of zt. This result implies that choosing the beamforming vectors 
to lie in the null spaces of the eavesdropper's channel vector, i.e., having | 2^m=i z mW m \ 2 = z^w = 0, 
is asymptotically optimal in the high-SNR regime. In this case, the eavesdropper cannot receive any data 
from the relays, and secrecy is automatically guarantied. No secrecy coding is needed at the relays. This 
asymptotic optimality can be seen from the following discussion. Assume that we impose the constraint 
z^w = 0. Now, the optimization problem (under the assumption No — 1) becomes 

2 \ 



max log 

w' w<Pt 
zt w =0 



V 



1 + 



2 4m— \ Z m W m 



J 



max log ( 1 + 

w<Pt 
zt w =0 



M 



max log I 1 

ztw=0 



m=l 
M 

m=l 



'P 



T 



log(Pr) + max log I J — 
wtw<i \ V Pr 

ztw=0 



M 



log (Pr) + l°g ( max 

wtw<l 
^ ztw=0 



M 



52 h * 

m=l 

2 N 



log(Pr) + log (max |h'N/>| 



m=l 
2\ 



(18) 

(19) 
(20) 

(21) 
(22) 



such that z^ = and H^ll 2 = 1- Above in (fl~9l) , we have defined w = w/y/Pr for which the constraint 
becomes w^w < 1. The approximation in (|2TT) is due to the fact that becomes negligible for large 
P T . Hence, null space beamforming provides the same asymptotic performance as in (flTT) and is optimal 
in the high-SNR regime. 
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Furthermore, the optimal null space beamforming vector can be obtained explicitly. Due to the null 
space constraint, we can write w = H^v, where denotes the projection matrix onto the null space of 
7) . Specifically, the columns of Hj; are orthonormal vectors that form the basis of the null space of z^. In 
our case, is an M x (M — 1) matrix. The power constraint wW = v^H^H^v = vV < Pt- Then, 
the optimization problem can be recast as 



max log I 1 + 

wt W<Py 



M 



5> 



:Wr. 



log ( 1 + max (w'hh'w} 

wtw<Pr 



log I 1 + max MH^hh'H^v 

\ vtv<P T 

log (l + P T A max (H^hhtH^)) 
log (l + P T tfHjHj f h) . 



(23) 

(24) 

(25) 
(26) 



Therefore, the optimum null space beamforming vector w is 



(27) 



where ci is a constant that is introduced to satisfy the power constraint. 

2) Low-SNR Regime: In the low SNR regime, in which both P s , Pt — > 0, we can see that 



v R\ . \9m 

lim — = mm 

P a ^0 P s m=X,...,M N„ 



and 



,. R s 
lim — 

Pa^O P T 



A max (hh t - zz f 



(28) 
(29) 



Thus, in the low SNR regime, the direction of the optimal beamforming vector approaches that of the 
eigenvector that corresponds to the largest eigenvalue of hh^ — zz^. A similar result is shown in a multiple- 
antenna setting in [[24]] . 



B. Optimal Beamforming under Individual Power Constraints 

In a multiuser network such as the relay system we study in this paper, it is practically more relevant to 
consider individual power constraints as wireless nodes generally operate under such limitations. Motivated 
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by this, we now impose \w m \ 2 < p m Vm or equivalently |w| 2 < p where | • | 2 denotes the element-wise 
norm-square operation and p is a column vector that contains the components {p m }. In what follows, the 
problem of interest will be again be the maximization of the secrecy rate or equivalently the maximization 
of the term inside logarithm function of R s d8]) but now under individual power constraints: 

Nn + I \^ A/ /l H) I 2 

max : ° ',^ =1 '„ (30) 



iwp< P n + 1 YZ 



2 



Nq + w'hbJ w 

= max — 1 — i — . (31) 

|w| 2 <p JMq + W T ZZ T W 

1) Semidefinite Relaxation (SDR) Approach: We first consider a semidefinite programming method 
similar to that in |fT9l . Using the definition X = ww\ we can rewrite the optimization problem in (|3D) 

as 

A^ + tr(hhtX) 

max T7 7 T^TV 

x A^o + tr(zztX) 

s.t diag(X) < p ^ 



rank X = 1, and X >- 



or equivalently as 



max t 
x.,t 



s.t tr(X(hh f - tzz f )) > N {t - 1), 
diag(K) < p, 

rank X = 1, and X >- 



(33) 



where tr(-) represents the trace of a matrix, diag(X) denotes the vector whose components are the diagonal 
elements of X, and X y means that X is a symmetric positive semi-definite matrix. The optimization 
problem in (l33l) is not convex and may not be easily solved. Let us now ignore the rank constraint in (l33l) . 
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That is, using a semidefinite relaxation (SDR), we aim to solve the following optimization problem: 

max t 
x,t 

s.t tr(X(hh f - tzz f )) > N {t - 1), ( 34 ) 

and diag(X) < p, and X >z 0. 

If the matrix X opt obtained by solving the optimization problem in (|34l) happens to be rank one, then 
its principal component will be the optimal solution to the original problem. Note that the optimization 
problem in (|34l) is quasiconvex. In fact, for any value of t, the feasible set in (|34|) is convex. Let i max be 
the maximum value of t obtained by solving the optimization problem (|34|) . If, for any given t, the convex 
feasibility problem 

find X 

such that tr(X(hh f - izz f )) > N (t - 1), (35) 

and diag(X.) < p, and X >z 

is feasible, then we have t max > t. Conversely, if the convex feasibility optimization problem d35l) is 
not feasible, then we conclude t max < t. Therefore, we can check whether the optimal value t mgx of the 
quasiconvex optimization problem in (|34l) is smaller than or greater than a given value t by solving the 
convex feasibility problem d35l) . If the convex feasibility problem (1351) is feasible then we know t max > t. 
If the convex feasibility problem (1351) is infeasible, then we know that t max < t. Based on this observation, 
we can use a simple bisection algorithm to solve the quasiconvex optimization problem (|34l) by solving 
a convex feasibility problem (|35l) at each step. We assume that the problem is feasible, and start with an 
interval [I, u) known to contain the optimal value t max . We then solve the convex feasibility problem at 
its midpoint t — (I +u)/2 to determine whether the optimal value is larger or smaller than t. We update 
the interval accordingly to obtain a new interval. That is, if t is feasible, then we set / = t, otherwise, we 
choose u — t and solve the convex feasibility problem again. This procedure is repeated until the width of 
the interval is smaller than the given threshold. Note that the technique of using bisection search to solve 
the SDP feasibility problem is also given in 11271 . Once the maximum feasible value for t max is obtained, 
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one can solve 



min trfX 
x 



s.t tr(X(hh^ — t max zz 



■ f )) > N (t 



max 



1) 



(36) 



and diag(X.) < p, and X y 



to get the solution X opi . (l36l) is a convex problem which can be solved efficiently using interior-point 
based methods. 

To solve the convex feasibility problem, one can use the well-studied interior-point based methods as 
well. We use the well-developed interior point method based package SeDuMi II341 . which produces a 
feasibility certificate if the problem is feasible, and its popular interface Yalmip [|35l . In semidefinite 
relaxation, the solution may not be rank one in general. Interestingly, in our extensive simulation results, 
we have never encountered a case where the solution X opt to the SDP problem has a rank higher than one. 
In fact, there is always a rank one optimal solution for our problem as will be explained later. Therefore, 
we can obtain our optimal beamforming vector from the principal component of the optimal solution X opi . 

2) Second-order Cone Program (SOCP) Approach: The reason that the SDR method is optimal for 
the above problem is that we can reformulate it as a second order cone problem E51 ll26ll by ignoring 
the phase in which we optimize w directly rather than performing the optimization over X = ww^. This 
provides us with another way of solving the optimization. The optimization problem (l30l) is equivalent to 



max t 

w,t 



(37) 



N + IbJwl 2 
N + Iztwl 2 



> t 



(38) 



and |w| 2 < p. 



Note that (l38l) can be written as 




(39) 
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where the equality on the right hand side of (|39l) follows from the definition of the magnitude- square of a 
vector. The equivalence of (f38|) and d39l) can easily be seen by rearranging the terms in (|39l) . In the above 
formulation, we have implicitly assumed that t > 1. Note that this assumption does not lead to loss of 
generality as we are interested in cases in which ^^|~T~p > 1- If this ratio is less than 1, the secrecy 
rate, as discussed before, is zero. 

Observe that an arbitrary phase rotation can be added to the beamforming vector without affecting the 
constraint in (|38l) . Thus, hW can be chosen to be real without loss of generality. We can take the square 
root of both sides of (|39l) . The constraint becomes a second-order cone constraint, which is convex. The 
optimization problem now becomes 



max t 



s.t W-h f w> 



/ 



(40) 



and |w| 2 < p. 



As described in the SDR approach, the optimal solution of (|4Q|) can be obtained by repeatedly checking the 
feasibility and using a bisection search over t with the aid of interior point methods for second order cone 
program. Again, we use SeduMi together with Yalmip in our simulations. Once the maximum feasible 
value t max is obtained, we can then solve the following second order cone problem (SOCP) to obtain the 
optimal beamforming vector: 



mm 



s.t 



w 



-bJw > 



z'w 



(41) 



and |w| < p. 



Thus, we can get the secrecy rate R s ,ind for the second-hop relay beamforming system under individual 
power constraints employing the above two numerical optimization methods. Then, combined with the 
first-hop source relay link rate Ri, secrecy rate of the decode and forward collaborative relay beamforming 
system becomes R dofMd = min(i?i, R s , ind )- 
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3) Simplified Suboptimal Design: As shown above, the design of the beamformer under individual 
relay power constraints requires an iterative procedure in which, at each step, a convex feasibility problem 
is solved. We now propose a suboptimal beamforming vector that can be obtained without significant 
computational complexity. 

We choose a simplified beamformer as w sim = 9w opt where w opt is given by (TT3T) with 1 1 w opt | | 2 = Pt = 
pt where pi is the individual power constraint for the i th relay, and we choose 



where w opt) k and pu are the /cth entries of w opt and p respectively, and we choose k as 

k = arg max - — opt ' 1 ■ (43) 
l<i<M pi 

Substituting this beamformer w sim into ©, we get the achievable suboptimal rate under individual power 
constraints. 

III. Amplify-and-Forward Relaying 

Another common relaying scheme in practice is amplify-and-forward relaying. In this scenario, the 
received signal at the m th relay R m is directly multiplied by l m w m without decoding, and forwarded to 
D. The relay output can be written as 
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1 



\Wopt,k\/y/Pk 



(42) 



x 




(44) 



The scaling factor, 



1 



(45) 



in 



^\g m \ 2 P s + Nj 
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is used to ensure E^av^l 2 ] = |u> m | 2 . The received signals at the destination D and eavesdropper E are 
the superposition of the messages sent by the relays. These received signals are expressed, respectively, as 



M 



Vd = ^2 h mW m l m (g m x s + rj m ) + no, and (46) 

m=l 
M 

y e = ^2 z m w m l m (g m x s + rim) + n x . (47) 



m=l 



Now, it is easy to compute the received SNR at D and E as 



IT h a I w \ 2 P 



EM j 1 2 p 

Em=l \Zm\ 2 l 2 rn \ W m\ 2 N m + N 

The secrecy rate is now given by 



(49) 



R s = I(x s ; y d ) - I(x s ; y e ) (50) 
= log(l + r d )-log(l + r c ) (51) 

I l m \ W m\ N m + N 



I I M j 12 p | Af _ : 2 70 I , 2 ^ ? ^ ' 

Em=lkm| 2 4k m | 2 iV m + iVo ' 



(52) 



Again, we maximize this term by optimizing {w m } jointly with the aid of perfect CSI. It is obvious that 
we only have to maximize the term inside the logarithm function. Let us define 

(53) 
(54) 
(55) 
(56) 



hg 


= [h\g\h,. 


••■> ^*m9*m^m] T , 


h z 


= [zlgth,.. 




D h 


= Diagd/i! 


| 2 Z 2 A^i,...,|/im| 2 4^), and 


D z 


= Diag(|zi| 


2 / 2 iY" l5 \zm\ 2 ImN m ). 
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Then, the received SNR at the destination and eavesdropper can be reformulated, respectively, as 

_ P s wth g h g t w _ P,tr(h g h g W) 
d wtD h w + iV ir(D h wwt) + JV ' 

P s wth z h z f w P s tr(h z h z 
e ~~ wtD z w + N ~ tr(D z wwt) + N 

With these notations, we can write the objective function of the optimization problem as 



, . r 1 I P S wthghgt w 

1 + 1 d wtD h w+AT 



1 + T P 1 i P a wth z h z T w 
- 1 wtD z w+AT 



(59) 



w+DhW + N + PsW+hghgW _ w f D z w + N 



wtD z w + N + P s wth z h z f w wtD h w + N 
N + tr((D h + P.hghgVwt) N + tr(D z ww 



x (60) 

(61) 



AT + tr((D z + P s h z h z t )wwt) iy + tr(D h wwt)' 
If we denote t, = ^o+M(D h+ p a h g h t) ww t) = jy +tr(D,wwt) and use h similar SDR h as 

described in the DF case, we can express the optimization problem as 

max tit 2 

X,tl,t2 

s.t tr(X(D z -f 2 D h )) > N (h - 1) 

(62) 

tr (X (D h + P.hghgt - h (D z + P s h z h z f ))) > No{h - 1) 

and diag(X) < p, (and/or tr(K) < Pt) and X y 0. 

Notice that this formulation is applied to both total relay power constraint and individual relay power 
constraint which are represented by ir(X) < P T and diag(X) < p, respectively. When there is only total 
power constraint, we can easily compute the maximum values of t\ and t 2 separately since now we have 
Rayleigh quotient problems. These maximum values are 

ti, u = \ max (r> h + + P.hghgt, D z + ^1 + P 8 hja/j , (63) 
t 2 , u = \ max (d z + ^1, D h + . (64) 
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When there are individual power constraints imposed on the relays, we can use the bisection algorithm 
similarly as in the DF case to get the maximum values t\^ u and £2,1," 3 f° r ^i an d ^2 by repeatedly solving 
the following two feasibility problems: 



find X 

s.t tr (X (D h + Pshgh^ - h (D z + P s h z h z f ))) > JV (*i - 1) (65) 
and diag(X.) < p, and X y 0, 



and 



find X 

s.t tr (X (D z - t 2 D h )) > jVo(t 2 -l) ( 66 ) 

and diag(K) < p, ararf X ^ 0. 

Note that for both total and individual power constraints, the maximum values of £1 and t 2 are obtained 
separately above, and these values are in general attained by different X = ww^. Now, the following 
strategy can be used to obtain achievable secrecy rates. For those X values that correspond to t\^ u and £ 1)M 
(i.e., the maximum t x values under individual and total power constraints, respectively), we can compute the 
corresponding t 2 = ^°+tr(Dhwwtj an d denote them as £ 2 ,m and i 2 ,« for individual and total power constraints, 
respectively. Then, log(£i i j jU t 2 and log(ti u i 2 /) will serve as our amplify-and-forward achievable rates 
for individual and total power constraints, respectively. With the achievable rates, we propose the following 
algorithm to iteratively search over t x and t 2 to get the optimal t l o and £ 2 ,o that maximize the product t\t 2 

Subscripts i in ti,i, u and t2,i,u are used to denote that these are the maximum values in the presence of individual power constraints. 
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by checking following feasibility problem. 

find X y 

s.t tr (X (D z - t 2 D h )) > N (t 2 - 1) 

fr (X (D h + P s h B hJ - U (D z + P.h.h.t))) > iV (ti - 1) ( 6V ) 
and tr(X) < Pr if there is total power constraint, 
or diag(X.) < p if there is individual power constraint. 

A. Proposed Algorithm 

Define the resolution At = or At = t -±jf- for some large iV for total and individual power constraints, 
respectively. 

1) Initialize t ljQ = t\ tU , t 2 , = h,i when total power constraint is imposed, and £i )0 = ti jijU , t 2 , = h,i,i 
when individual power constraint is imposed. Initialize the iteration index i = N. 

2) Set ti = zAt. If tit 2iU < ti ;0 t 2; o (total power constraint) or tit 2jijU < ti )0 t 2)0 (individual power 
constraint), then go to Step (3). Otherwise, 

a) Let t 2 = 1, ° t2, ° . Check the feasibility problem (|67l) . If it is infeasible, % = i — 1 go to step (2). 
If it is feasible, use the bisection algorithm in (1671) with t\ to get the maximum possible values 
of t 2 and denote this maximum as t 2 rn . The initial interval in the above bisection algorithm 
can be chosen as l ' 1 '^ 2 '" , t 2)U ] or [ * 1 ^ 2 -° ,t 2t i }U ] depending on the power constraints. 

b) Update ti = t\, t 2 o = t 2 m , % — i — 1. Go back to step (2). 
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3) Solve the following problem to get the optimal X 

min tr(X) 
x 

s.t tr (X (D z - t 2)0 D h )) > N (t 2<o - 1) 
tr (X (D h + P fl h g h g t - t lj0 (D z + PshJaJ))) > N {t liO - 1) 

(68) 

X y and 

tr(X) < P T if there is total power constraint, 
diag(X.) < p if there is individual power constraint. 

B. Discussion of the Algorithm 

Our algorithm is a two-dimensional search over all possible pairs (ti, £ 2 )» which can produce the greatest 
product t±t 2 , whose logarithm will be the global maximum value of the secrecy rate. In the following, we 
will illustrate how our algorithm works for individual power constraints. Similar discussion applies to the 
total power constraint case as well. The algorithm initiates with the achievable pair *2,i,0' m which 

ti i U is the maximum feasible value for ti. Thus, all t\ values in our search lie in [0, tiiu]. We chose the 
resolution parameter N to equally pick N points in this interval. We then use a brute force strategy to 
check each point iteratively starting from ti jijU down to 0. In each iteration, the feasibility problem (1671) 
is quasi-convex. Thus, we can use the bisection search over t 2 to get the greatest value of t 2 . Note that 
our initial bisection interval for t 2 is [ fl, °* 2 '° , t 2yiyU ] where t 2 ,i,« is the maximum feasible value for t 2 , and 
tl; °' 2, ° is chosen so that the optimal t 2 we find at the end of the bisection search will produce a product 
t\t 2 that is greater than our currently saved optimal ti^t^o- With this approach, after each iteration, if a 
t 2 value is found, the new optimal t\ Q t 2 o will be greater than the previous one. Note that our iteration's 
stop criterion is t\t 2 ^ u < t^ t 2y0 . This means that further decrease in the value of t\ will not produce a 
product tit 2 that is greater than our current t\ Q t 2 o . Thus, the value t\ Q t 2 o at the end of this algorithm will 
be the global maximum since we have already checked all possible pairs t\,t 2 that are candidates for the 
optimal value. 
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IV. Robust Beamforming Design 

All of the beamforming methods discussed heretofore rely on the assumption that the exact knowledge 
of the channel state information is available for design. However, when the exact CSI is unavailable, the 
performance of these beamforming techniques may degrade severely. Motivated by this, the problem of 
robust beamforming design is addressed in ll32l and [|33l . The robust beamforming for MISO secrecy 
communications was studied in ||2~9ll where the duality between the cognitive radio MISO channel and 
secrecy MISO channel is exploited to transform the robust design of the transmission strategy over the 
secrecy channel into a robust cognitive radio beamforming design problem. 

We additionally remark that, beside the assumption of perfect channel state information, our previous 
analysis is applicable only when the relays are fully synchronized at the symbol level. When the time 
synchronization between the relays is poor, the signal replicas passed through different relays will arrive 
to the destination node with different delays. This will result in inter-symbol-interference (IS I). To combat 
such ISI, the authors of ll30l view an asynchronous flat-fading relay network as an artificial multipath 
channel (where each channel path corresponds to one particular relay), and use the orthogonal frequency 
division multiplexing (OFDM) scheme at the source and destination nodes to deal with this artificial 
multipath channel. In OTTl . a filter-and-forward protocol has been introduced for frequency selective relay 
networks, and several related network beamforming techniques have been developed. In these techniques, 
the relays deploy finite impulse response (FIR) filters to compensate for the effect of source-to-relay and 
relay-to-destination channels. Since the relay synchronization problem is out of the scope of this paper, 
we will mainly focus on combatting the effect of imperfect channel state information in the following 
discussion. 

Systems robust against channel mismatches can be obtained by two approaches. In most of robust 
beamforming methods, the perturbation is modeled as a deterministic one with bounded norm which lead 
to a worst cast optimization. The other approach applied to the case in which the CSI error is unbounded 
is the statistical approach which provides the robustness in the form of confidence level measured by 
probability. 

Let us consider the DF case. We define H = hh^ and Z = zz-t as the channel estimators, and H = H — H 

20 



and Z = Z — Z as the estimation errors. First, consider the worst case optimization. In the worst case 
assumption, H and Z are bounded in their Frobenius norm as ||H|| < e#, ||Z|| < e^, where en,£z are 
assumed to be upper bounds of the channel uncertainty. Based on the result of |f32]| , the robust counterpart 
of previously discussed SDR-based optimization problem can be written as 

max t 

s.t ir(X((H - e H l) - i(Z + e z l)) > N (t - 1), ( 69 ) 

and diag{X) < p, and X >z 0. 

Note that the total power constraint tr(X) < P T can be added into the formulation or substituted for the 
individual power constraint in (|69| ). This problem can be solved the same way as discussed before. 

However, the worst-case approach requires the norms to be bounded, which is usually not satisfied in 
practice. Also, this approach is too pessimistic since the probability of the worst-case may be extremely low. 
Hence, statistical approach is a good alternative in certain scenarios. In our case, we require the probability 
of the non-outage for secrecy transmission is greater than the predefined threshold e by imposing 

Pr ( *o + M(H + H)X) >f \ =pr A / X(II + fi _ t{± + ^ > N (t _ > £ (7Q) 
\iV + tr((Z + Z)X) " y V V "~ '))- 

Now, the optimization problem under imperfect CSI can be expressed as 

max t 

s.t Pr (tr (X(H + H - i(Z + Z)) > N (t - 1))) > e, ( 71 ) 

arac/ diag{X) < p (or tr(X) < P^), and X ^ 0. 

If relays are under individual power constraints, we use diag(X) < p. Otherwise, for the case of total 
power constraint, we use tr(X) < P T . We can also impose both constraints in the optimization. 

Note that the distribution of the components of the error matrices H and Z depend on the channel 
estimation technique and distribution of the channel coefficients. In order to simplify the analysis and 
provide an analytically and numerically tractable approach, we assume that the components of the Hermitian 
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channel estimation error matrices H and Z are independent, zero-mean, circularly symmetric, complex 
Gaussian random variables with variances o 2 ^ and cr|. Such an assumption is also used in 11331 . Now, we 
can rearrange the probability in the constraint as 

Pr (tr ((H-tZ + U- tZ)X\ > (t - l)N j\ . (72) 

Let us define y = tr f(H — tZ + H — tZ)Xj. For given X, H, and Z, we know from the results of 
11331 that y is a Gaussian distributed random variable with mean \i = tr ^(H — tZ)xj and variance 
a 2 = (er^ + t 2 cr|) tr(XX' 1 '). Then, the non-outage probability can be written as 

Pr(y >(t- 1)N ) = — exp (-^Ll^\ (73) 

J{t-i)N v2na y V la y J 



or equivalently as, 

(^-l)iVo-^ errl( _ 2£ + 1) (75) 

V2(Jy 

Note that s should be close to one for good performance. Thus, both — 2e + 1 and should be 



negative valued. Note further that we have tr(XX"l") = ||X|| 2 , and hence a y = yo"^- + ^ 2cr |||X||. Then, 
this constraint can be written as 

|x|| < ('-W-J' . (76) 

^(al + tVperf-^^e + l) 
As a result, the optimization problem becomes 

max t 
x,t 

s.t iixii < c-w-/. 



'2(4 + fV|)err 1 (-2,r+l) 
and diag(X.) < p(or tr(X) < Py), and X >z 0. 
Using the same bisection search, we can solve this optimization numerically. 



(77) 
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V. Numerical Results 

We assume that {g m }, {h m }, {z m } are complex, circularly symmetric Gaussian random variables with 
zero mean and variances a 2 , a\, and a 2 z respectively. We first provide numerical results for decode-and- 
forward beamforming schemes. In our numerical results, we focus on the performance of second-hop 
secrecy rate since the main emphasis of this paper is on the design of the beamforming system in the 
second-hop. Moreover, each figure is plotted for fixed realizations of the Gaussian channel coefficients. 
Hence, the secrecy rates in the plots are instantaneous secrecy rates. 

In Figures [2] and [3l we plot the second-hop secrecy rate, which is the maximum secrecy rate that our 
collaborative relay beamforming system can support under both total and individual relay power constraints. 
For the case of individual relay power constraints, we assume that the relays have the same power budgets: 
Pi = -^p. Specifically, in Fig. [2l we have cr^ = 3, a z = 1, Nq = 1 and M = 5. In this case, the legitimate 
user has a stronger channel. In Fig. [3l the only changes are cx^ = 1 and a z = 2, which imply that the 
eavesdropper has a stronger channel. Our CRB system can achieve secure transmission even when the 
eavesdropper has more favorable channel conditions. As can be seen from the figures, the highest secrecy 
rate is achieved, as expected, under a total transmit power constraint. On the other hand, we observe that 
only a relatively small rate loss is experienced under individual relay power constraints. Moreover, we note 
that our two different optimization approaches give nearly the same result. It also can be seen that under 
individual power constraint, the simple suboptimal method suffers a constant loss as compared to SDR or 
SOCP based optimal value. 

In Fig. HI we fix the relay total transmitting power as Pt = lOdB, and vary the number of collaborative 
relays. Other parameters are the same as those used in Fig. [3j We can see that increasing M, increases 
the secrecy rate under both total and individual power constraints. We also observe that in some cases, 
increasing M can degrade the performance when our simplified suboptimal beamformer is used. 

In Fig. [51 we plot the secrecy rate for amplify-and-forward collaborative relay beamforming system for 
both individual and total power constraints. We also provide the result of suboptimal achievable secrecy 
rate for comparison. The fixed parameters are a g = 10, ah = 2,a z = 2, and M = 10. Since the AF secrecy 
rates depend on both the source and relay powers, the rate curves are plotted as a function of Pt/P s - As 
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Fig. 3. DF Second-hop secrecy rate vs. the total relay transmit power Pt for different cases. Eavesdropper has a stronger channel. 
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Fig. 4. DF second-hop secrecy rate vs. number of relays for different cases. 



before, we assume that the relays have equal powers in the case in which individual power constraints are 
imposed, i.e., pi = P T /M. It is immediately seen from the figure that the achievable rates for both total 
and individual power constraints are very close to the corresponding optimal ones. Thus, the achievable 
beamforming scheme is a good alternative in the amplify-and-forward relaying case due to the fact that 
it has much less computational burden. Moreover, we interestingly observe that imposing individual relay 
power constraints leads to only small losses in the secrecy rates with respect to the case in which we have 
total relay power constraints. 

In Fig. [6l we plot the maximum second hop secrecy rate of decode-and-forward that we can achieve for 
different power P T and non-outage probability e values. In this figure, we fix M = 5. h and z are randomly 
picked from Rayleigh fading with o~ h = 1 and o- z = 2, and we assume that estimation errors are inversely 
proportional to P T . More specifically, in our simulation, we have a 2 ^ = 0.1/ P T and <r| = 0.2/ P T . We also 
assume the relays are operating under equal individual power constraints, i.e., p^ = It is immediately 
observed in Fig. [6] that smaller rates are supported under higher non-outage probability requirements. In 



25 



12 



o 

.Q 

E 
>■ 




4 







5 



10 



15 



20 



PJP 

T s 



Fig. 5. AF secrecy rate vs. P T /Ps- cr g = 10, a h = 2, a z = 2, M = 10. 

particular, this figure illustrates that our formulation and the proposed optimization framework can be used 
to determine how much secrecy rate can be supported at what percentage of the time. For instance, at 
Pt = 20dB, we see that approximately 7.4 bits/symbol secrecy rate can be attained 70 percent of the time 
(i.e., e = 0.7) while supported secrecy rate drops to about 6.2 bits/symbol when e = 0.95. 



In this paper, collaborative beamforming for both DF and AF relaying is studied under secrecy constraints. 
Optimal beamforming designs that maximize secrecy rates are provided under both total and individual relay 
power constraints. For DF with total power constraint, we have remarked that the optimal beamforming 
vector is the solution of a Rayleigh quotient problem. We have further identified the beamforming structure 
in the high- and low-SNR regimes. For DF with individual relay power constraints and AF with both total 
and individual relay power constraints, we have formulated the problem as a semidefinite programming 
problem and provided an optimization framework. We have also provided an alternative SOCP method to 
solve the DF relaying with individual power constraints. In addition, for DF relaying, we have described 



VI. Conclusion 
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Fig. 6. DF second secrecy rate vs. Pt under different e. 



the worst-case robust beamforming design when CSI is imperfect but bounded, and the statistical robust 
beamforming design based upon minimum non-outage probability criterion. Finally, we have provided 
numerical results to illustrate the performance of beamforming techniques under different assumptions, 
e.g., DF and AF relaying, total and individual relay power constraints, perfect and imperfect channel 
information. 
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